First-order phase transitions from poles in asymptotic 
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Although partition functions of finite-size systems are always analytic, and hence have no poles, 
they can be expressed in many cases as series containing terms with poles. Here we show that such 
poles can be related to linear branches of the entropy, expressed in the thermodynamic limit as a 
function of the energy per particle. We also show that these poles can be used to determine whether 
the entropy is nonconcave or has linear parts, which is something that cannot be done with the 
sole knowledge of the thermodynamic free energy derived from the partition function. We discuss 
applications for equilibrium systems having first-order phase transitions. 
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It is known from the seminal work of Lee and Yang [T] 
that the zeros of partition functions, seen as functions of 
a complex variable, provide useful information about the 
analyticity properties of the free energy in the thermo- 
dynamic limit, and thus about the appearance of phase 
transitions in that limit. Their main results, extended by 
Fisher [2 to the canonical partition function, 

microstatcs uj 

of an A^-particle system with Hamiltonian H, show that 
the complex zeros of Z(jS) lie away from the real axis for 
all N < oo, and that, in the presence of a phase transi- 
tion, the zeros of Z(j3) get infinitesimally close to the real 
axis in the limit N — > oo. The real value j3 c at which the 
zeros "pinch" the real axis in this limit corresponds to the 
critical inverse temperature at which the thermodynamic 
free energy function, defined as 

<p(J3) = lim -^lnZ(/3), (2) 

is non-analytic [3]. Moreover, the angle at which the 
locus of zeros crosses the real axis determines the order 
of the phase transition corresponding to (3 C [3]. These 
properties of the zeros of Z(j3) have been widely used for 
studying phase transitions in equilibrium systems (see [5] 
for a review), as well as, more recently, nonequilibrium 
systems in steady states [6]. 

We study in this paper another component of Z(f3), 
which we refer to as the set of "poles" of Z(/3), and re- 
late it to the thermodynamic properties of the system 
represented by H. We use the term "pole" with quotes 
because Z{j3) has of course no singular points, since it is 
an analytic function of j3 for N finite. The poles that we 
have in mind actually appear in the terms of a special 
asymptotic expansion of Z(j3) to be defined below. The 
main result that we prove here is that these poles are re- 
lated to linear branches of the microcanonical entropy of 



the system represented by H, and can therefore be used 
to study two physical phenomena associated with these 
branches, namely, first-order phase transitions and phase 
separation in the canonical ensemble [7]. 

The knowledge of these poles also solves an outstand- 
ing problem in the field of long-range interacting systems 
[8], which is to determine whether a system undergoing 
a first-order phase transition in the canonical ensemble 
has an entropy in the microcanonical ensemble which is 
nonconcave or is concave but has a linear branch. These 
two types of entropy are known to lead to the same non- 
differentiable free energy function ip(f3) [7], so they can- 
not be distinguished from the point of view of this func- 
tion alone. Our results show, however, that they can 
be distinguished using information about the "poles" of 
Z(f3). As an illustration of these results, we compute the 
linear entropy of a simple model of DNA denaturation 
from its partition function. At the end, we also com- 
ment on the applicability of our results for calculating 
large deviation functions characterizing the fluctuations 
of nonequilibrium systems. 

The problem that we are concerned with is to calculate 
the microcanonical entropy function 

s(u) = lim — lnO(u) (3) 

from the density of states which is obtained from 

Z(f3) via the inverse Laplace transform formula 

i />r+ioo 

n{u) = — Z([3)e? Nu dp, (4) 

i7Tl J, — ioo 

where r is an arbitrary real number located inside the 
region of convergence of Z(/3) [S]. Given that s(u) is a 
thermodynamic-limit quantity, it is often assumed that 
its calculation via Eq. {£} requires not the exact knowl- 
edge of Z(f3), but only of the asymptotic behavior of 
Z((3) as N — > oo, expressed, according to Eq. as 
Z(/3) » e~ N<fi W with sub-exponential corrections in N. 
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By substituting this asymptotic result in the inverse La- 
place transform, and by performing a saddlepoint ap- 
proximation of the complex integral, one indeed finds 

s(u)=wi{Pv,-<p(P)}. (5) 

The problem with this result, which is nothing but a 
Legendre transform written in a technical form, is that 
it does not always hold [7]. To see why, consider the 
following two "mock" partition functions: 

_ -AT/3 

Z 1 (P)=e N P + e- N f l , Z 2 {P) = - j- . (6) 

It is easily verified that each of these partition functions 
is analytic, despite the appearance of the 1//3 term in 
the second, and that both partition functions lead to the 
same free energy function (f(/3) = However, the 

densities of states obtained from Z\(JS) and Z 2 {(3) and 
their corresponding entropies must be different, since the 
partition functions are themselves different. A simple 
calculation of the inverse Laplace transform shows that 
this is the case. The density of states obtained from 
Zi((3) is fli(u) = S(u+1) + S(u—1), whereas the density 
of states obtained from Z 2 (f3) is such that £l 2 (u) = 1 for 
u 6 [—1,1] and is otherwise. As a result, the entropy 
s\{u) obtained from VL\(u) is finite and zero only for u — 
±1, whereas the entropy S2(u) associated with tt 2 (u) is 
finite and zero for all u £ [—1, 1]. Hence the entropies 
calculated from Z\(ji) and Z 2 (j3) are different, but both 
partition functions lead to the same free energy (p(/3). 

The problem illustrated by this example is well docu- 
mented in large deviation theory [7J, and has been dis- 
cussed recently in the context of long-range interaction 
systems, such as self-gravitating particles and unscreened 
plasmas, as these often have nonconcave entropies [5]. 
The problem is that entropy functions that have the same 
concave envelope, as in the example above, have the same 
free energy </?(/?), and cannot, as mentioned before, be 
distinguished from the point of view of </?(/3) alone. This 
means concretely that the knowledge of the asymptotic 
behavior Z{fi) ss e ~ Ni f{P) j s no t sufficient in general to 
compute s(u); extra information is required to faithfully 
obtain s(u), especially if one suspects that s(u) is non- 
concave or has a linear branch. 

The two partition functions shown in Eq. Q give a 
hint as to what kind of extra information is required. By 
recasting each of these partition functions in the form 

Z{p)=a{p)e N ^+b{P)e- N P, (7) 

we see that what distinguishes Zi(j3) from Z 2 (f3) is the 
presence of poles in the coefficients a(/3) and b((3) of 
Z 2 {j3). As we show next, it is the presence of these poles 
in the series representation of Z 2 (f3) that is responsible 
for the linear behavior of s 2 (u) seen for u e [—1, 1]. This 
applies to any partition function, in the sense that poles 



in asymptotic expansions of Z(f3) are generally associated 
with linear branches of s(u). 

To demonstrate this claim, we go back to the formula 
of the inverse Laplace transform shown in Eq. Q, and 
assume, as an extension of Eq. Q, that Z{j3) admits an 
asymptotic expansion of the form 

Z{fi)=Y t c iW)e-"' pm - ( 8 ) 
j 

This expansion can always be obtained for ID systems 
by expanding Z(/3), for example, in the eigenbasis of the 
transfer matrix associated with H [TO]. For systems of 
higher dimensions, there is not necessarily a transfer ma- 
trix to work with, and for these, Eq. ^ should presently 
be considered as an ansatz rather than a derived result. 
This point will be discussed in more detail in a subse- 
quent paper [llj . 

Here we shall work on the assumption that Eq. ^ is 
given, and that the functions ipj obtained are concave 
and smooth functions of (3 that do not depend on N. 
Moreover, we shall assume that the coefficients Cj((3) are 
sub-exponential in N, and may have poles in the complex 
/3-plane [T2]. These assumptions are verified for some 
models of interest [10], including the one studied at the 
end of this paper. 

It should be mentioned that expansions similar to 
Eq. ^ have been considered before in studies of Yang- 
Lee zeros, first-order phase transitions, and metastability 
(see, e.g., [T3] )• However, to our knowledge, none have 
considered the possibility that the coefficients Cj (/3) may 
have poles in /3. To see how these poles relate to the 
properties of Q(u) and in turn s(u), we insert the expan- 
sion Q in the formula of the inverse Laplace equation, 
and proceed to evaluate the complex integral by going 
through the following steps: 

1 . Distribute the integral of the inverse Laplace trans- 
form inside the sum of the partition function to obtain 

= £ ^ / z mu ~ vm] d/3. (9) 

This is permitted provided that all the integrals inside 
the sum converge. With this in mind, we should choose 
r in such a way that each integral in Eq. ([9]) is convergent. 
In particular, we cannot put r on any poles of Cj(f3). 

2. Approximate each of the integrals labeled by j in 
Eq. ([9]) to exponential order in N using the saddlepoint or 
steepest descent approximation [14] . This requires that 
we deform the vertical integration contour sitting at r, 
which is often called the Bromwich contour, to another 
equivalent steepest-descent contour that passes through 
the saddlepoint of the exponent 

$ i G8,«) = /3«-p j ( j 9), (10) 

in such a way that Im$j(/3, it) is constant. Assuming 
that <fij((3) is differentiable and concave, the saddlepoint 
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is given by the unique solution of f'AP) 



Hence- 



forth, we denote this saddlepoint by /?* , and the steepest- 
descent contour passing through this point by Dj. 

3. Two situations will arise from the previous step, 
depending on whether or not the deformation of the 
Bromwich contour to the steepest-descent contour neces- 
sitates that we cross poles of Cj(f3). On the one hand, if 
no such poles need to be crossed, then the integral on the 
Bromwich contour is equivalent to the integral evaluated 
on Dj. On the other hand, if the deformation requires 
that we cross any poles of Cj (/?) , then 
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Res ' 



(ii) 



where Res is the sum of the residues of the poles that 
were crossed when transforming the Bromwich contour 
B into the steepest-descent contour Dj (see Fig. [IJ. 

4. Insert the result obtained in pTj ) into the sum of 
Eq. ([9]), and cancel any terms that have the same mag- 
nitude but opposite sign. At this point, we expect many 
residue terms to cancel. What remains can be put in the 
form 

^) = Efi / ^(^e^^d/S + ^Res^V 

(12) 

where Res(/3^) denotes the residue of Cj(f3) e w *3'(0> u ) for 
the pole (3* e . 

5. Approximate the integrals on the steepest-descent 
contours Dj by their saddlepoints, i.e., 



2?ri 



Cj (j3)e N *^dj3 « e^W'") (13) 



D , 



with sub-exponential corrections in N [13] ■ The term 
Cj(p) does not contribute to the approximation because 
it is assumed to be sub-exponential with N. 

6. Evaluate the residue terms. If we assume, for sim- 
plicity, that the Cj(f3)'s have only simple poles, then the 
residues are approximately given by 



Res(/3 



(14) 



with sub-exponential corrections in N, where <jji = ±1. 
Note that this approximation must be performed only for 
those poles P* t that were crossed in Step 3 and do not 
get cancelled in Step 4 above. The remarkable feature 
of these poles is that they give rise to terms that are 
exponential in N similar to the saddlepoints. 

7. The result of Steps 5 and 6 is the following approx- 
imation for the density of states: 



fi(u) 



+ > tr jt e 1(P *' 



(15) 




FIG. 1. Deformation of the Bromwich contour B sitting at r 
to a steepest-descent contour D crossing the saddlepoint /3*. 
A residue arises if the deformation crosses a pole /3 X . 



As a last step, we further approximate this expression by 
taking the largest exponential term (Laplace approxima- 
tion). In order to express this final step in a convenient 
form, let us define Bj to be the set containing the sad- 
dlepoint (3* and the poles /3, x „ that remain after Step 4. 
Then, by taking the largest term in that expression, we 
obtain 



s(u) = sup sup {j5u — ifj(/3)}. 



(16) 



Note that the sign variable <jji does not appear in the 
above result because the dominant term of £l(u) is nec- 
essarily positive. Moreover, the saddlepoint or pole se- 
lected from the maximization over Bj is necessarily real, 
for otherwise s(u) would not be a real function. 

The representation of the entropy s(u) shown in 
Eq. (16) is the main result of this paper. The concav- 



ity properties of s(u) are deduced from this equation by 
analyzing the maximization over the set Bj. Three cases 
must be distinguished: 

Case 1 : If, for an interval U of values for u, the maxi- 
mization over Bj picks up a pole, then s(u) will be pro- 
portional to u over hi (since poles of Cj(/3) do not depend 
on u). In this case, s(u) will thus have a linear branch. 

Case 2: If, for u G hi, the maximization over Bj does 
not pick up a pole, but picks up instead a saddlepoint /3* 
which is constant as a function of u, then s(u) will also 
have a linear branch over hi. 

Case 3: If, for u £ hi, the maximization over Bj picks 
up neither a pole (Case 1) nor a saddlepoint j3* which 
is constant in u (Case 2), then s(u) will vary in a non- 
linear way with u. In this case, s(u) may be concave or 
nonconcave, but has no linear branch over hi. 

It can be checked that the linear entropy s^iu) men- 
tioned earlier arises from Case 1 above. The pole in the 
corresponding partition function Z 2 {/3) is at /3 = and 
gives rise to the flat part of s(u) with slope 0. Case 3 
applies, on the other hand, to Z\{fi), as the latter has no 
pole in its series representation. For examples of entropy 
calculations involving constant saddlepoints, see [T5" l ll6 | . 

To provide an explicit illustration of Case 1 involving 
poles, we now calculate the entropy of a simple model of 
DNA denaturation due to Kittel [T7] . The partition func- 
tion of this model can be written in the thermodynamic 
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limit as 



i 



-N(Pe-lnG) 



(17) 



where e is the energy associated with one bond of a DNA 
chain consisting of N bonds, G is a degeneracy factor as- 
sociated with this energy, and (3 C = (T 1 In G is the critical 
inverse temperature at which the thermodynamic free en- 
ergy <p(/3), derived from Eq. ([2|, is nondifferentiable [18]. 
From the Yang-Lee theory point of view, j3 c is also the 
accumulation point of the zeros of Z(f3). 

To obtain s(u) for this model, we first note that Z(f3) 
in Eq. ^TtJ) has the form of Eq. pH) with (pi(fi) = and 
f2{P) = (3e — InG, and that each of the two terms com- 
posing Z(/3) has a pole at (i c . Next we choose r > (3 C , 
and follow the calculation steps described above. For 
u G [0, e), we find that the Bromwich integral involv- 
ing the first term of Z(f3) gives rise to a residue pro- 
portional to e N P" u because of the pole at /3 C , whereas 
the Bromwich integral of the second term vanishes [19] . 
Therefore, s(u) = (3 c u for u £ [0, e). For all other values 
of u, the Bromwich integrals either cancel one another or 
vanish, and so we find s[u) — -co outside [0, e). This 
agrees with the entropy that one would obtain by com- 
binatorial means. Moreover, the fact that the finite part 
of s(u) is linear with slope (3 C confirms the fact that <£>(/?) 
is nondifferentiable at f3 c [7]. 

This calculation of s(u), although simple, provides an 
illustration of what should be observed in more realis- 
tic equilibrium systems with first-order transitions, espe- 
cially systems involving short-range interactions, such as 
nearest-neighbor spin systems or screened Coulomb sys- 
tems [5]. For these, it is known that s(u) is in general a 
concave function in the thermodynamic limit |20j . This 
implies that, if <p(/3) is nondifferentiable, then s(u) will 
have in general one or more linear branches [7J, which 
are likely to be associated, according to our results, with 
poles in some expansion of Z{(3). 

It should be obvious, to conclude, that our results can 
be applied to partition functions other than the canonical 
one to calculate the entropy as a function of macrostates 
other than the energy per particle (e.g., magnetization 
or particle density). Our results can also be generalized, 
following the theory of large deviations [7J, to calculate 
entropy functions describing the fluctuations of observ- 
ables of nonequilibrium systems in driven steady states. 
In this context, one must replace Z(f3) by the generat- 
ing function of the observable considered. Poles in series 
representations of generating functions have been consid- 
ered in the context of nonequilibrium systems (see, e.g., 
[T51 [H]), and are known to be associated with exten- 
sions of the Gallavotti-Cohen symmetry of nonequilib- 
rium fluctuations [22] . We expect such poles to also play 
a role in nonequilibrium first-order phase transitions, as 
these are generally characterized by nonequilibrium en- 
tropy functions (viz., rate functions) having linear or non- 



convex branches. 
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